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Prediction  of  the  frequency-dependent  dielectric  function  of  thin  films  poses  computational  chal¬ 
lenges,  and  at  the  same  time  experimental  characterization  by  spectroscopic  ellipsometry  remains 
difficult  to  interpret  because  of  changes  in  stoichiometry  and  surface  morphology,  temperature, 
thickness  of  the  him,  or  substrate.  In  this  work,  we  report  calculations  for  titanium  nitride  (TiN),  a 
promising  material  for  plasmonic  applications  because  of  less  loss  and  other  practical  advantages 
compared  to  noble  metals.  We  investigated  structural,  electronic,  and  optical  properties  of  stoichio¬ 
metric  bulk  TiN,  as  well  as  of  the  TiN(100),  TiN(llO),  and  TiN(lll)  outermost  surfaces.  Density 
functional  theory  (DFT)  and  many-body  GW  methods  (Green’s  (G)  function-based  approximation 
with  screened  Coulomb  interaction  (W))  were  used,  ranging  from  GoWo,  GWo  to  partially  self- 
consistent  sc-GWq,  as  well  as  the  GW-BSE  (Bethe-Salpeter  equation)  and  time-dependent  DFT 
(TDDFT)  methods  for  prediction  of  the  optical  properties.  Structural  parameters  and  the  band 
structure  for  bulk  TiN  were  shown  to  be  consistent  with  previous  work.  Calculated  dielectric  func¬ 
tions,  plasma  frequencies,  reflectivity,  and  the  electron  energy  loss  spectrum  demonstrated  consis¬ 
tency  with  experiment  at  the  GW0-BSE  level.  Deviations  from  experimental  data  are  expected  due 
to  varying  experimental  conditions.  Comparison  of  our  results  to  spectroscopic  ellipsometry  data 
for  realistic  nanostructures  has  shown  that  although  TDDFT  may  provide  a  computationally  feasi¬ 
ble  level  of  theory  in  evaluation  of  the  dielectric  function,  application  is  subject  to  validation  with 
GW-BSE  calculations.  ©  2015  A1P  Publishing  LLC.  [http://dx.doi.org/10. 1063/1.4935813] 


I.  INTRODUCTION 

Titanium  nitride  gained  much  interest  in  the  last  decades 
because  of  unique  properties, 1-3  such  as  high  hardness,  high 
wear  and  corrosion  resistance,  which  are  useful  for  protec¬ 
tive  and  decorative  coatings,4'5  compatibility  with  comple¬ 
mentary  metal-oxide-semiconductor  (CMOS)  processes, 
e.g.,  used  as  diffusion  barriers  and  gate  material,6  or  for  solar 
cell  application.7  Recently,  TiN  thin  films  were  suggested  as 
alternative  plasmonic  materials,8-12  exhibiting  less  loss  and 
providing  other  practical  advantages  compared  to  noble  met¬ 
als.  For  example,  comparing  TiN  and  gold  disk-shaped  nano¬ 
particles,  it  was  shown  that  the  TiN  nanodisks  are  efficient 
local  heat  sources  and  outperform  gold  in  the  biological 
transparency  window.11  Furthermore,  fabrication  of  epitaxial 
superlattices  with  thin  film  TiN  as  a  plasmonic  component 
for  optical  hyperbolic  metamaterials  demonstrated  a  higher 
photonic  density  of  states  than  the  gold-  or  silver-based 
respective  nanomaterial  structures.12  In  order  to  assess  the 
bounds  of  performance  of  such  materials,  knowledge  of  the 
material  properties  is  required  in  solving  Maxwell’s  equa¬ 
tions,  e.g.,  with  the  finite-difference  time-domain  method.1’ 
In  this  case,  a  best  fit  of  an  experimental  frequency- 
dependent  dielectric  function  to  an  analytic  expression  is 
used  in  the  simulated  spectral  region. 


a,Author  to  whom  correspondence  should  be  addressed.  Electronic  mail: 
ruth.pachter@us.af.mil 


The  frequency-dependent  dielectric  function  of  TiN  thin 
films  was  characterized  by  spectroscopic  ellipsometry  (SE) 
in  previous  work,8'14-21  analyzed,  for  example,  with  a 

2  f  2 

Drude-Lorentz  model  e  =  - -  'f - b  'vL,  2 — 1 2°J  -r  1 

00  io2-n  D<o  ‘—‘J=  l 

Soo  is  a  background  constant  due  to  interband  transitions, 
larger  than  unity  due  to  higher-energy  transitions  not  repre¬ 
sented  by  the  i  Lorentz  oscillators  considered.  In  the  Drude 
term,  mpu  is  the  unscreened  plasma  frequency,  correlated 
with  the  free  electron  density  and  the  metallic  character,  and 
To  is  a  damping  factor  due  to  the  scattering  of  electrons, 
determined  by  the  inverse  of  the  mean  time  between  colli¬ 
sions,  i.e.,  free  electron  relaxation  time.  The  /  Lorentz  oscil¬ 
lators  have  interband  excitation  m0j,  oscillator  strength  f, 
and  damping  factors  r;.  Various  parameters  were  varied  to 
fit  the  SE  data  for  TiN,  e.g.,  the  number  of  Lorentz  oscilla¬ 
tors,22  where  models  with  two  or  four  Lorentz  oscillators16 
were  adopted.14  Indeed,  parameters  that  change  the  response 
of  the  film  are  generally  not  well  defined,  such  as  effects  of 
deposition  conditions,  including  stoichiometry,  surface  mor¬ 
phology  and  crystallinity,  temperature,  bias  voltage,  or  sub¬ 
strate  and  thickness.  For  instance,  the  unscreened  plasma 
frequency  was  shown  to  decrease  substantially  depending  on 
the  deposition  temperature  or  nitrogen  flow  rates,  while 
effects  of  substrate  could  be  important  for  films  thinner  than 
60nm.14  Below  5nm,  the  dielectric  functions  change  drasti¬ 
cally  with  decreasing  film  thickness.20  In  the  work  of  Naik 
et  a/.,12  5nm  thick  layers  were  used.  In  another  example,  it 
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was  demonstrated  that  stoichiometric  TiN(OOl)  layers  grown 
on  MgO(OOl)  are  fully  relaxed  at  the  growth  temperature, 
but  changes  are  noted  for  sub-stoichiometric  structures, 
where  lattice  constants  decrease  for  TiNx  (x  as  low  as  0.7) 
grown  on  Mg(001)  at  700°C,23  which  in  turn  will  affect  the 
optical  properties. 

At  the  same  time,  although  density  functional  theory 
(DFT)  calculations  of  structural  and  electronic  properties  of 
TiN  bulk  and  surface  structures  were  previously  per¬ 
formed.2  1  analysis  of  the  optical  properties  at  an  ab  initio 
level  for  this  material  was  not  addressed.  The  GW  approxi¬ 
mation  within  the  framework  of  many-body  perturbation 
theory36,37  emerged  as  useful  in  prediction  of  optical  proper¬ 
ties  using  the  BSE  method  to  include  excitonic  effects.38^0 
GW-BSE  (Bethe-Salpeter  equation)  is  computationally  in¬ 
tensive,  but  linear  response  time-dependent  DFT  (TDDFT) 
could  be  inadequate  in  predicting  optical  spectra  of  solid 
materials  and  depends  on  the  exchange-correlation  func¬ 
tional  applied.41  However,  notwithstanding  the  experimental 
variations  expected  and  the  computational  limitations,  it  is  of 
interest  in  the  first  stage  to  obtain  a  reference  frequency- 
dependent  dielectric  function  theoretically  by  accurate  first 
principles  calculations  for  stoichiometric  TiN. 

In  this  work,  we  report  calculations  of  the  frequency- 
dependent  real  and  imaginary  parts  of  the  dielectric  function 
(e  =  eW  +  ie1'2'1)  for  stoichiometric  bulk  TiN,  as  well  as  for 
surfaces  that  mimic  the  outermost  layer,  in  comparison  to  SE 
data.  The  goal  is  to  assess  prediction  of  the  complex  refrac¬ 
tive  index  in  the  low-energy  region  resulting  from  intraband 
transitions  (lower  than  ca.  2  eV  (Ref.  14)),  and  from  inter¬ 
band  transitions  at  higher  energy.  For  this  purpose,  we  per¬ 
formed  a  systematic  study  of  structural,  electronic,  and 
optical  properties  for  stoichiometric  bulk  TiN,  and  TiN(lOO), 
TiN(llO),  TiN(lll)  surfaces,  applying  DFT/TDDFT  and 
GW  methods.  Dielectric  functions  were  calculated  using 
GW-BSE  and  TDDFT  with  generalized  gradient  approxima¬ 
tion  (GGA)  and  range-separated  exchange-correlation  func¬ 
tionals.  We  found  that  surface  geometries  and  band 
structures  were  consistent  with  previous  work.  In  bench¬ 
marking,  the  GWq-BSE  method  was  shown  to  be  appropriate 
for  the  larger  systems.  In  addition,  comparison  of  the  GWq- 
BSE  and  TDDFT  dielectric  functions  has  shown  that  the 
functional  is  subject  to  validation  with  GW-BSE  calcula¬ 
tions.  Analysis  of  the  GWq-BSE  results  in  comparison  to  SE 
data  demonstrated  reasonable  agreement  with  experiment, 
where  variations  due  to  experimental  conditions  have  to  be 
accounted  for.  This  investigation  will  lead  to  further  charac¬ 
terization  of  TiN  thin  films  and  provide  a  basis  for  the  inter¬ 
pretation  of  measured  optical  properties. 

II.  METHODS 

Calculations  were  performed  using  the  Vienna  ab  initio 
simulation  package  VASP  5.3  (Refs.  42  and  43)  within 
the  framework  of  DFT.  The  projector  augmented-wave 
(PAW)  method  was  applied4’  for  describing  the  electron-ion 
interaction,  and  the  Perdew-Burke-Emzerhof  (PBE)  GGA 
functional44  was  used,  unless  otherwise  indicated.  The  Kohn- 
Sham  equations43,46  were  solved  using  a  plane-wave  basis 


set,  with  an  energy  cutoff  of  400  eV.  For  the  surface  calcula¬ 
tions,  a  supercell  with  a  vacuum  region  in  the  three- 
dimensional  cell  was  used.  Convergence  for  an  adequate 
vacuum  size  was  tested,  and  values  of  10-16  A  were  found  to 
be  large  enough  to  ensure  that  the  periodically  repeated 
images  do  not  interact.  Atomic  positions  were  relaxed  until 
forces  on  each  atom  were  smaller  than  0.01  eV/A. 
Monkhorst-Pack  grid47  sampling  of  the  Brillouin  zone  (BZ) 
with  6x6x6  k  points  for  the  bulk  TiN  and  6x6x1  k 
points  for  surfaces  was  applied.  A  series  of  independent  con¬ 
vergence  tests  that  corroborate  these  parameters  to  provide 
the  total  energies  within  a  few  meV  was  performed.  The  cut¬ 
off  energy  was  kept  fixed  when  performing  convergence  test¬ 
ing  of  the  number  of  k  points,  size  of  vacuum,  number  of 
surface  layers,  and  relaxation  algorithm  used.  Mixing  and 
smearing  parameters  were  chosen  to  obtain  improved  agree¬ 
ment  with  experiment. 

TiN  crystalizes  in  the  NaCl  lattice  structure  within  the 
Fm3m  space  group.  Three  facets  of  the  TiN  surface  were 
modeled,  namely,  TiN(100),  TiN(llO),  and  TiN(l  1 1),  as  well 
as  bulk  TiN.  TiN(100)  and  TiN(llO)  surfaces  were  modeled 
by  ten  TiN  layers  with  Ti  and  N  atoms  in  the  same  plane 
(total  of  twenty  atoms).  Ti  and  N  atoms  in  the  TiN(lll) 
surface  are  not  in  the  same  plane,  but  at  a  perpendicular  dis¬ 
tance  of  1.22  A  from  each  other.  This  leads  to  two  possibil¬ 
ities  for  the  TiN(lll)  surface,  namely,  Ti-  or  N-terminated 
structures.  Twenty  one  layers  of  TiN(lll)  with  ten  and 
eleven  layers  of  either  Ti-  or  N-termination  were  used.  The 
(1  x  1)  cell  of  surfaces  is  shown  in  Fig.  1.  Optimization  of 
TiN(100),  TiN (11 0),  and  Ti-  and  N-terminated  TiN(lll) 
surfaces  was  carried  out  using  PBE.  The  equilibrium  lattice 
constant  of  bulk  TiN,  calculated  by  minimization  of  the  total 
energy  of  the  system  with  respect  to  the  total  volume  of  the 
unit  cell,  resulted  in  an  optimized  lattice  constant  of  4.25  A, 
in  agreement  with  the  experimental  values  of  4.24  A 
(Ref.  48)  and  4.25  A  (Ref.  49)  for  bulk  TiN,  thus  validating 
the  functional  used.  Our  result  is  also  consistent  with  previous 
work  using  GGA,  with  values  of  4.24-4.26  A.26,27,30-33  The 
lattice  constant  of  bulk  TiN  calculated  with  LDA  (local  den¬ 
sity  approximation),  of  4.175  A,  is  underestimated  compared 
to  the  GGA  values.26  Note  that  these  values  decrease  for  x 
smaller  than  1  in  TiNx.23 

Band  structures  were  calculated  by  DFT  and  GW  meth¬ 
ods.  Details  on  the  implementation  of  the  GW  method  within 
the  PAW  framework  are  given  by  Shishkin  and  Kresse  in 
Refs.  50  and  51.  GW-BSE39  and  TDDFT52  were  used  for 
calculation  of  the  frequency-dependent  dielectric  function. 
The  large  computational  cost  associated  with  GW  calcula¬ 
tions  required  the  use  of  fewer  surface  layers,  and  the  system 
size  was  reduced  to  5  layers  for  the  TiN(100),  TiN(llO) 
surfaces,  and  11  layers  for  TiN(  1 1 1)  in  this  case.  The  itera¬ 
tive  approach  for  GW  calculations  was  carried  out  as 
follows.  A  non-self-consistent  G0W0  evaluation  of  the  quasi¬ 
particle  self-energy  was  performed  on  top  of  the  computa¬ 
tionally  less  demanding  DFT  GGA  calculation.  Next,  only  G 
was  iterated  and  W  kept  fixed  (GW0),  generally  shown  to 
provide  better  agreement  with  experimental  data,51  and  is 
relatively  computationally  less  expensive.  To  perform  self- 
consistent  GW0  (sc-GWq)  calculations,  non-diagonal 
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FIG.  1.  Top  and  side  views  of  TiN(lOO) 
(a),  TiN(llO)  (b),  and  N-terminated 
TiN(l  1 1)  (c)  surfaces. 


components  of  the  self-energy  are  included  in  the  orbital  ba¬ 
sis,  only  iterating  G  and  keeping  W  fixed.  sc-GW0  is  the 
most  computationally  intensive  of  the  three  approximations 
we  applied,  and  the  /.-point  mesh  was  reduced  to  4  x  4  x  1 . 
GW  calculations  that  fully  update  both  Green’s  function  G 
and  the  screened  potential  W,  namely,  full  self-consistent 
GW  (sc-GW),  are  computationally  very  intensive  and  their 
utility  is  still  controversial,41'51  while  GW0  or  sc-GW0  are 
commonly  used.53,54  Calculation  of  the  dielectric  function  is 
carried  out  as  a  four  step  process.  In  the  first  step,  a  ground 
state  DFT  calculation  is  performed  for  occupied  orbitals,  and 
in  the  second  stage,  the  wavefunction  from  the  previous  step 
is  used  to  perform  a  ground  state  calculation  with  a  large 
number  of  unoccupied  empty  bands  and  the  dielectric  matrix 
is  calculated.  In  the  third  step,  GW  calculations  (GoW0,  GW0, 
or  sc-GWq)  are  performed  by  using  the  wavefunction  and 
dielectric  matrix  from  the  previous  step,  kept  fixed  throughout 
the  various  steps.  At  this  stage,  the  band  structure  is  also  cal¬ 
culated.  Band  structures  were  calculated  using  both  non-self- 
consistent  and  self-consistent  evaluation  of  the  quasiparticle 
self-energy.  In  the  last  step,  the  BSE  is  solved’9  to  obtain  the 
dielectric  function. 


TABLE  I.  wpu  values  (eV)  for  bulk  TiN  and  TiN  surfaces  (eV). 


(fipu 

GWo-BSE 

TDDFT/PBE 

TDDFT/HSE 

Bulk  TiN 

7.62 

6.12 

6.72 

TiN(100) 

4.76 

3.60 

3.97 

TiN(llO) 

3.50 

3.76 

4.16 

TiN(l  1 1)-Ti  terminated 

3.34 

2.05 

2.29 

TiN(l  1 1)-N  terminated 

2.84 

1.39 

1.56 

The  imaginary  part  of  the  frequency-dependent  dielec¬ 
tric  function  (s  =  £f  1  -1  +  itP1 )  in  the  PAW  method  (following 
the  notation  of  Gajdos  et  al ,55)  is  given  by  the  Cartesian 

tensor  s^(m)  =  lim^o  ^  J2c.v.k2w^(wkeck  -  Gk  -  w) 
X  (uck+eaq\uvk)  (uck+eal,q\uYk}*,  where  0)  is  in  units  of  energy; 
Q  the  volume  of  the  primitive  cell;  e7  unit  vectors  for  the 
three  Cartesian  directions;  the  indices  c  and  v  refer  to  con¬ 
duction  and  valence  states,  respectively;  Mrk  are  the  cell- 

(2) 

periodic  orbitals  at  k.  In  the  calculation  of  (to),  k  is 
restricted  to  the  irreducible  wedge  of  the  BZ.  The  real  part  of 
the  dielectric  tensor  is  obtained  from  the  Kramers-Kroning 

transformation,  given  by  s^j  (w)  =  1  +  =PJ0  J_co2+iridco, 
where  P  denotes  the  principal  value  of  the  integral  and  the 
small  complex  shift ;/  (CSHIFT)  is  used  for  smoothing  of  the 
spectrum,  which  was  set  to  0.02  eV  in  all  cases.  Intraband 


FIG.  2.  Plasma  frequency  as  a  function  of  the  number  of  k  points  for  bulk 
TiN  and  TiN(  100),  TiN(l  10)  and  N-  and  Ti-terminated  TiN(l  1 1)  surfaces. 
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TABLE  II.  Percentage  surface  relaxation  of  TiN(lOO),  TiN(llO),  and 
TiN(lll)  surfaces  with  respect  to  bulk-terminated  distances  (A);  T  and  D 
for  TiN(lOO)  and  TiN(llO)  surfaces,  as  explained  in  the  text.  Negative  and 
positive  signs  indicate  contraction  and  expansion  between  the  surfaces, 
respectively. 


TiN(100) 

TiN(llO) 

TiN(lll) 

T 

D 

T 

D 

N-terminated  Ti-terminated 

di,2 

18.1 

-1.10 

15.02 

-11.04 

-36.63 

-8.39 

^2,3 

0.59 

0.45 

-8.76 

4.70 

25.92 

0.28 

^3,4 

-1.58 

0.77 

-2.88 

0.17 

-7.15 

2.78 

d4,5 

2.22 

0.46 

0.92 

0.97 

1.28 

2.77 

ds.6 

0.69 

1.37 

-0.89 

1.57 

-2.11 

^6.7 

-0.41 

0.39 

^7,8 

0.23 

0.85 

d8,9 

0.82 

-0.29 

^9,10 

-1.41 

1.07 

transitions  were  calculated  for  the  metallic  system,56'57 

where  the  real  part  is  expressed  as  e^1"  *™(co)  =  1  — 
with  the  square  of  the  intraband  plasma  frequency  given  by 

»’«/;  =  ^E„,k2gk%^(ea%^)  (e/J^r)-  The  frequency- 
dependent  optical  constants  n  and  k  were  derived  from  the 

calculated  dielectric  function  by  n—  W^(g(''  +  \J g^^  +  g®2) 

and  k=\J |(— gT)  4-  yj e(i)2_|_£(2^  as  n+ik=  \J gW  +  ig®, 

and  the  electron  energy  loss  spectrum  (EELS)  by 
gt2)/(g(')2  -)-  g(“)7) .  The  maximum  frequency  for  the  dense 
part  of  frequency  grid  (OMEGAMAX)  was  set  to  lOOeV, 
and  the  number  of  frequency  grid  points  (NEDOS)  to  a  large 
value  of  10000  for  a  dense  set  of  data  points.  The  number  of 
bands  (NBAND)  was  set  to  64  for  TiN  bulk  and  surfaces, 
providing  a  significant  number  of  empty  bands.  The  band 
convergence  was  tested  by  increasing  the  number  to  256,  but 
no  significant  differences  were  noted;  conventionally,  the 


FIG.  3.  Band  structures  and  LDOS  of  bulk  TiN  and  TiN(100),  TiN(l  10),  N-  and  Ti-terminated  TiN(l  11)  surfaces. 
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number  of  bands  is  set  to  half  the  number  of  electrons  in  the 
system,  plus  the  total  number  of  atoms,  plus  a  few  additional 
empty  bands.  Earlier  studies  suggested  that  reasonably  con¬ 
verged  results  could  be  obtained  by  choosing  bands  that  are 
1.5  times  the  number  of  electrons.1'8 

The  /.--point  meshes  were  tested,  specifically  4x4x4, 
5x5x5,  and  6x6x6,  and  as  small  differences  were 
noted,  a  6x6x6  /.-point  mesh  for  GoW0  and  GW0  and 
4x4x4  for  sc-GW0  calculations  were  selected.  A  larger 
mesh  size  greatly  increases  memory  and  computational 
requirements  when  solving  the  BSE  and  is  not  practical. 
Other  studies  using  similar  methods  also  found  that  such 
/-point  meshes  provide  spectra  in  good  agreement  with  the 
experimental  data.59  The  number  of  empty  bands  and  the 
size  of  the  /-point  mesh  affect  performance  of  the  calcula¬ 
tion  by  making  it  very  demanding  computationally,  and  our 
choice  assured  accurate  prediction.  We  found  that  a  very 
dense  /-point  mesh  is  essential  to  obtain  a  well-converged 
plasma  frequency.  Calculations  of  the  plasma  frequency  with 
4x4x4  (4x4x1)  to  40  x  40  x  40  (40  x  40  x  1)  /-point 
meshes  for  the  bulk  (surface)  were  performed.  In  most  cases, 
we  found  the  plasma  frequency  to  reasonably  converge 
around  24  x  24  x  24  (24  x  24  x  1)  Z  points  (see  Table  I  and 
Fig.  2).  However,  such  a  dense  mesh  is  not  feasible  for  GW- 
BSE  calculations  because  of  the  large  memory  requirements. 
Therefore,  a  sparser  mesh  was  used  to  obtain  the  interband 
part  of  the  dielectric  function.  For  comparison,  TDDFT  cal¬ 
culations  of  the  dielectric  function  using  the  PBE  or  the 
range-separated  Heyd-Scuseria-Ernzerhof  (HSE06,  range- 
separation  parameter  0.11  bohr  1  (Refs.  60  and  61))  func¬ 
tionals  were  performed.  Ground  state  calculations  were 


performed  in  this  case  using  large  /-point  meshes  of 
12x  12x  12  and  12  x  12  x  1  for  the  bulk  and  surface  sys¬ 
tems,  respectively,  followed  by  TDDFT  calculations. 

III.  RESULTS  AND  DISCUSSION 
A.  Surface  morphology 

TiN  surfaces  are  known  to  undergo  rippling,  character¬ 
ized  by  deviations  in  the  interlayer  separation  from  that  of 
the  bulk  material,  and  by  the  distance  between  neighboring 
surface  planes.  In  bulk-terminated  surfaces  prior  to  optimiza¬ 
tion,  the  perpendicular  distance  is  zero,  as  the  atoms  lie  in  a 
plane.  During  structural  relaxation,  non-zero  di  i+1  distances 
(T)  and  average  displacement  deviation  of  the  planes  (D) 
result.  T  and  D  values  are  summarized  in  Table  II,  reported 
as  a  percentage  change,  where  positive  and  negative  signs 
denote  expansion  or  contraction  between  the  surfaces, 
respectively.  Bulk-terminated  distances  between  two  surface 
planes  in  TiN(100)  and  TiN(llO)  are  2.12A  and  1.50 A, 
respectively.  As  seen  in  Table  II,  surface  formation  for 
TiN(100)  does  not  perturb  the  equilibrium  distances  signifi¬ 
cantly.  Variation  in  T  for  the  top  layer  compared  to  the 
bulk-terminated  surface  (18.1%  or  0.18  A)  is  consistent  with 
the  17.9%  change  that  was  previously  calculated.26'32 
Displacement  T  values  are  less  than  2%  for  d2i3  and  d3  4, 
lower  than  previously  calculated,  where  a  smaller  number  of 
layers  was  however  used.26'32  Surface  relaxation  quantified 
by  D  (see  Table  II)  shows  a  small  contraction  of  ca.  — 1%, 
while  the  rest  of  the  TiN  (100)  surface  layers  expand  slightly 
for  compensation.  Earlier  calculations26  reported  a  small 
contraction  of  about  —2%,  and  an  expansion  of  about  5% 


s 


(2) 


FIG.  4.  Calculated  dielectric  functions 
e  as  a  function  of  the  frequency  (eV). 


5 

Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 


1 95302-6 


Mehmood  etal. 


J.  Appl.  Phys.  118,  195302  (2015) 


was  also  reported/2  The  rippling  effect  of  the  TiN(llO)  sur¬ 
face  is  similar  to  that  of  TiN(lOO),  with  T  between  the  first 
and  second  layer  of  0.15  A  (or  15%),  decreasing  for  subse¬ 
quent  interlayer  distances,  comparable  to  previously 


calculated  values  of  0.146-0.163  A.26’31,32  A  relatively  large 
contraction  in  the  distance  between  the  first  and  second  sur¬ 
face  planes  (D  in  Table  II)  was  calculated,  of  —11%,  result¬ 
ing  in  an  expansion  of  about  5%  between  the  second  and 


e(1)  -  GW0-BSE 

-12  -10  -8  -6-4-2  0  2 


s'”  -  GW0-BSE  8«  -  GW0-BSE 


FIG.  5.  Dielectric  functions  s  calculated  at  various  levels  of  theory  in  comparison  to  GW0-BSE  for  TiN  bulk,  TiN(100)  and  TiN(llO)  surfaces  ((a)-(c)), 
respectively. 
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third  layers,  consistent  with  prior  work.26,31,32  Earlier  calcu¬ 
lations  also  showed  that  a  converged  value  for  the  relaxation 
is  found  for  a  range  of  9-1 1  atomic  layers.31 

The  N-terminated  TiN(lll)  surface  demonstrated  a 
large  contraction  between  the  first  and  second  layers, 


accompanied  by  a  smaller  expansion  for  the  second  and  third 
layers,  while  corresponding  changes  in  the  Ti-terminated 
TiN(lll)  surface  were  much  smaller,  as  was  similarly 
reported.  ’  ’  Note  that  the  surface  thickness  in  our  mate¬ 
rial  models  ranges  from  ~1.4nm  for  TiN(llO)  to  ~2.5  nm 


eV 


FIG.  6.  Calculated  dielectric  functions  for  bulk  TiN  and  TiN(lOO),  TiN(l  10),  and  N-  and  Ti-terminated  TiN(l  1 1)  surfaces,  as  a  function  of  frequency  (eV),  in 
comparison  to  SE  data  (a,14  b,20  c8)  for  (a)  GW0-BSE,  (b)  TDDFT-HSE,  and  (c)  TDDFT-PBE. 
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for  TiN(lll).  To  understand  effects  of  a  larger  number  of 
layers  on  surface  relaxation,  molecular  dynamics  simulations 
with  an  appropriate  force-field  are  required,  for  example,  as 
performed  by  Sangiovanni  et  al.  for  admolecule  transport  on 
TiN(OOl)  surfaces,62  to  be  addressed  in  future  work. 

B.  Electronic  band  structures 

Using  the  optimized  bulk  and  surface  TiN  structures, 
calculations  of  band  structures  were  performed,  as  shown  in 
Fig.  3.  Our  results  are  in  agreement  with  the  work  of  Yin 
et  al.,63  as  shown  Fig.  SI  in  the  supplementary  material.64 
Interband  transitions  for  bulk  TiN  at  ca.  1,  2.3,  3.0-3. 9,  and 
5.0-6.2eV  were  previously  calculated,14’16’28’65’66  and  gaps 
in  the  3.35-3.45  eV  range  were  reported. 67~69  The  local  den¬ 
sity  of  states  (LDOS)  of  bulk  TiN  and  the  top  layer  of 
TiN(lOO),  TiN(llO),  and  TiN(lll)  surface  facets  are  shown 
in  Fig.  3.  Three  main  regions  can  be  identified.  For  bulk 
TiN,  the  lowest  valence  band  comprises  N  5  states  and  some 
hybridized  Ti  d  states  around  — 8eV.  The  second  lowest 
sub-band  consists  of  strongly  hybridized  N  p  states  and  Ti 
d  states  representing  the  bonding  part  ranging  from  about 
—1  eV  to  the  Fermi  level  (OeV).  The  third  region,  the  conduc¬ 
tion  band,  consists  mainly  of  Ti  d  states  and  a  small  degree  of 
N  p  states  corresponding  to  the  antibonding  part.  When  surfa¬ 
ces  are  formed,  a  change  in  the  LDOS  is  noted.  A  small  shift 
of  the  center  of  the  band  (^0.5  eV)  towards  higher  energies 
for  the  Ti  atom  in  the  top  layer  was  calculated,  which 
decreases  when  moving  away  from  the  surface.  The  shift  is 
stronger  for  N  atoms  and  amounts  to  1  eV  for  TiN(100)  and 
ca.  2eV  for  TiN(lll)  surfaces.  Strongest  hybridization  is 
noted  for  the  N-terminated  TiN( 111)  surface. 

C.  Optical  properties 

Frequency-dependent  dielectric  functions  for  bulk  TiN, 
and  TiN(100),  TiN(llO)  and  TiN(lll)  surfaces  were  calcu¬ 
lated  using  the  GW-BSE  method  and  TDDFT  with  the  PBE 
and  FISE  functionals.  Comparison  of  the  results  is  shown  in 
Fig.  4  for  bulk  TiN  and  the  considered  surfaces.  We  com¬ 
pared  GWo-BSE  results  for  up  to  3  eV  to  those  using  other 
levels  of  theory  because  this  is  the  spectral  region  that  could 
be  relevant  to  plasmonic  applications.  For  analysis,  calcu¬ 
lated  GWq-BSE  fv ri  and  e!-2'’  (calculated  at  intervals  of 
0.125  eV),  in  comparison  to  GoW0-BSE,  sc-GW0-BSE  and 
TDDFT,  are  shown  in  Figs.  5(a)-5(c)  for  TiN  bulk  and 
TiN(100)  and  TiN(llO)  surfaces,  respectively  (results  for 
TiN(lll)  are  shown  in  Fig.  S264).  sc-GW0-BSE  dielectric 
functions  are  in  good  agreement  with  the  GW0-BSE  results, 
but  this  is  not  always  the  case  for  GoW0-BSE.  For  example, 
for  ^  for  the  TiN(100)  surface,  the  slope  in  the  linear 
regression  is  0.5  in  this  case,  but  changes  to  0.9-1  for  other 
levels  of  theory  (Fig.  5(b)).  For  bulk  TiN,  TDDFT  using 
PBE  and  FISE  functionals  provided  slopes  that  deviate  con¬ 
siderably  from  1,  for  example,  2. 6-2. 7  for  d2;  (Fig.  5(a)). 
GW0-BSE  eW  and  vs.  TDDFT  for  TiN(l  10)  demonstrate 
deviations,  with  correlation  coefficients  of  0.8-0. 9,  as  com¬ 
pared  to  values  close  to  1  for  bulk  TiN  and  TiN(100).  The 
agreement  of  TDDFT  with  GW0-BSE  for  the  TiN(l  11)  sur¬ 
face  is  even  worse  (see  Fig.  S2),  yet  still  to  be  considered 


because  of  the  computational  cost  of  GW-BSE  calculations. 
Thus,  calculated  dielectric  functions  at  the  GW0-BSE,  and 
TDDFT-HSE  and  TDDFT-PBE  levels  were  chosen  for  com¬ 
parison  with  SE  data. 

Next,  we  compare  GW0-BSE  dielectric  functions  for 
bulk  TiN  and  TiN(100),  TiN(llO),  and  TiN(lll)  surfaces  to 
experimental  data,  as  summarized  in  Fig.  6(a),  and  values 
obtained  using  TDDFT  with  HSE  and  PBE  functionals  in 
Figs.  6(b)  and  6(c),  respectively.  On  comparing  the  experi¬ 
mental  data  for  in  Ref.  8,  measured  at  300  °C,  to  the  data 
in  Ref.  14,  and  Ref.  20,  a  relative  consistency  is  demon¬ 
strated,  except  for  the  region  of  up  to  about  2eV,  where 
intraband  transitions  are  important,  and  where  the  plasma 
frequency  is  affected  by  deposition  conditions.  For  e1-1'1 ,  var¬ 
iations  in  the  experimental  data  are  larger,  and  so  are  differ¬ 
ences  in  the  theoretically  determined  values,  and  it  is 
therefore  difficult  to  compare  with  the  calculated  values 
directly.  Plasma  frequency  copu  values  are  listed  in  Table  I. 
Our  results  are  consistent  with  previous  calculations,70  and 
with  experimental  measurements,  however  dependent  on  the 
bias  voltage  and  temperature  during  synthesis.  Larger  volt¬ 
age  biases  (from  40  V  to  200  V)  or  temperatures  (up  to 
600  °C)  will  increase  mpu  from  ca.  5eV  up  to  7.5eV.14  SE 
measured  data  are  also  dependent  on  the  thickness,14’20  and 
the  stoichiometry  of  the  thin  film  with  nitrogen  vacancies  or 
atom  incorporation  in  interstitial  positions.71  TDDFT  results 
reproduce  less  well  experimental  data,  particularly  for  the 
N-terminated  TiN(lll)  surface.  Notably,  for  both  the 
and  the  N-terminated  TiN( 111)  surface  provided  poorer 


e'11  Experimental  Data 

-10  -8-6-4-2  0  2 


FIG.  7.  Calculated  dielectric  functions  s  vs.  experimental  data  from  Ref.  8. 
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FIG.  8.  Refractive  index  n  (a),  extinc¬ 
tion  coefficient  k  (b),  reflectivity  (c), 
and  EELS  spectrum  (d)  for  bulk  TiN  at 
the  sc-GW0-BSE  level  of  theory. 


general  agreement  with  experiment.  As  we  noted  in  examin¬ 
ing  the  surface  structure,  the  N-terminated  TiN(lll)  surface 
demonstrated  a  large  contraction  between  the  first  and  sec¬ 
ond  layers,  and  further  investigation  of  details  of  structural 
parameter  changes  of  the  surface  is  to  be  undertaken  in 
future  work. 

For  a  quantitative  comparison  with  GW0-BSE,  we  used 
the  data  in  Ref.  8,  as  shown  in  Fig.  7.  A  reasonable  agree¬ 
ment  with  experiment  is  demonstrated  for  bulk  TiN  for  efI ' 
and  £^  (correlation  coefficients  of  0.9-1).  However,  this  is 
generally  not  the  case  for  TiN  surfaces  as  expected,  because 
the  SE  data  we  compared  to  has  been  obtained  for  larger 
thicknesses  than  the  outermost  layers  modeled.  Validation  of 
the  theoretical  results  depends  on  careful  experimental  char¬ 
acterization  of  the  structure,  and  calculations  for  amorphous 
structures  will  be  carried  out  in  future  work.  Optical  con¬ 
stants,  namely,  refractive  index  n  and  extinction  coefficient  k 
for  bulk  TiN  calculated  at  the  sc-GW0-BSE  level,  are  shown 
in  Figs.  8(a)  and  8(b),  respectively.  The  theoretical  results 
overestimate  the  data  from  Palik,72  commonly  used  in  elec¬ 
tromagnetic  simulations  of  plasmonic  nanostructures,  partic¬ 
ularly  in  the  low-energy  spectral  region,  as  expected.  A 
better  agreement  is  noted  with  the  data  of  Pfliiger  et  a/.49 
The  reflectivity  spectrum  of  bulk  TiN  is  shown  Fig.  8(c), 
having  a  high  absolute  value  at  low  energies,  indicative  of  its 
metallic  behavior,  and  consistent  with  SE  measurements.16 
The  EELS  spectrum  in  Fig.  8(d)  is  also  consistent  with 
experiment  for  bulk  TiN.16 

IV.  CONCLUSION 

In  this  work,  we  have  undertaken  an  investigation  of 
electronic  and  optical  properties  of  stoichiometric  bulk  TiN, 
and  TiN(lOO),  TiN(llO),  TiN(lll)  surfaces,  benchmarking 
the  band  structure  and  comparing  the  dielectric  function  to 
SE  data  obtained  for  realistic  samples.  To  address  theoretical 
challenges  in  prediction  of  electronic  and  optical  properties 


of  solid  materials,  both  DFT  and  many-body  GW  methods 
were  applied,  ranging  from  GoW0  to  partially  self-consistent 
sc-GWq,  and  TDDFT  (using  the  PBE  and  HSE  functionals) 
and  GW-BSE  for  dielectric  function  predictions.  The  struc¬ 
tural  ripple  effect  in  TiN  was  quantified  by  optimization  of 
surface  structures,  and  the  optimized  structures  used  for  cal¬ 
culation  of  the  electronic  and  optical  properties.  The  band 
structure  we  calculated  demonstrated  consistency  with  previ¬ 
ous  work  for  bulk  TiN.  GW0-BSE  dielectric  functions,  ana¬ 
lyzed  in  comparison  to  SE  data  of  the  real  and  imaginary 

parts  of  the  dielectric  function,  demonstrated  qualitative 
agreement  with  experiment,  where  variations  due  to  experi¬ 
mental  sample  characteristics  have  to  be  considered. 
Consistency  with  experiment  was  also  shown  in  calculation 
of  the  plasma  frequency,  as  well  as  the  reflectivity  and  EELS 
spectra.  We  found  that  TDDFT  is  less  suitable  for  evaluation 
of  the  optical  constants  n  and  k,  yet  could  be  considered 
upon  benchmarking  with  higher  level  calculations.  A  calcu¬ 
lated  dielectric  function  for  materials  for  which  experimental 
characterization  is  complicated  can  be  used  in  electromag¬ 
netic  simulations  of  plasmonic  nanostructures  to  assess 
bounds  of  performance,  or  in  design  of  multi-layer  optical 
coatings  that  require  knowledge  of  the  frequency-dependent 
complex  refractive  index. 
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